AJ, ACCEPTED 2006/12/1 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



O 
O 

(N 

O 

Q 

00 
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ABSTRACT 

We present models of the shapes of four Kuiper belt objects (KBOs) and Jovian Trojan (624) 
Hektor as ellipsoidal figures of equilibrium and Roche binaries. Our simulations select those figures of 
equilibrium whose lightcurves best match the measured rotational data. The best fit shapes, combined 
with the knowledge of the spin period of the objects provide estimates of the bulk densities of these 
objects. We find that the lightcurves of KBOs (20000) Varuna and 2003 ELgi are well matched by 
Jacobi triaxial ellipsoid models with bulk densities 992j^g kg m -3 and 2551^J 5 kgm -3 , respectively. 
The lightcurves of (624) Hektor and KBO 2001 QG298 are well-described by Roche contact binary 
models with densities 2480 l 2 ^ 2 kg m~ 3 and 590t\^ 3 kgm~ 3 , respectively. The nature of 2000 GN m 
remains unclear: Roche binary and Jacobi ellipsoid fits to this KBO are equivalent, but predict 
different densities, ^2000 kg m~ 3 and ~650kgm -3 , respectively. Our density estimates suggest a 
trend of increasing density with size. 

Subject headings: Kuiper Belt objects — minor planets, asteroids — solar system: general 
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1. INTRODUCTION 

Most small bodies of the solar system appear unre- 
solved at the ~0.05 arcsecond peak angular resolution 
offered by current technology. As a consequence, infor- 
mation about the shapes and rotations of the small bod- 
ies must be inferred, principally from measurements of 
the time dependence of the scattered radiation. So-called 
"lightcurve inversion" techniques have been used for 
decades to study the rotational properties of main-bel t 
asteroids (jCellino et al.lfl989l : iKaasalainen et al.l 12002] ) . 
At their simplest, these involve using the peak-to-peak 
interval of the lightcurve to estimate the rotation period 
and the peak-to-peak brightness variation to estimate the 
"axis ratio" of bodies that are assumed to be triaxial in 
shape and in principal axis rotation about the minor axis. 
At their most complex, the inversion techniques can be 
used to solve for the full 3-dimensional shapes and rota- 
tion vec tors, utilizing observations from a range of aspect 
angle s (jKaasalainen fc Torppal 120011 : IKaasalainen et all 
2001). 

All lightcurve interpretations are subject to an am- 
biguity between variations caused by shape and varia- 
tions caused by non-unifo rm surface alb edo, as clearly 
expressed a century ago bv lRusseiH £l906). This ambigu- 
ity can be broken when simultaneous optical and thermal 
observations are available, as is the case for some of the 
larger asteroids in the main-belt. Numerous observations 
of this type have shown that, with rare exceptions, the 
albedos of the aster oids do not vary ove r their surfaces 
by a large amount (|Degewii et al.l ll979). This spatial 
uniformity could simply mean that the compositions are 
intrinsically uniform. Alternatively, real surface compo- 
sitional variations could exist but be smoothed-out by ef- 
ficient lateral transport of dust over the surfaces of small 
bodies. The most famous exception to this rule is pro- 
vided by Saturn's satellite Iapetus, which has a lightcurve 
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range of n early two ma gnitudes caused by surface albedo 
markings (Mi lIisllT977l ). This case is pathological, how- 
ever, in the sense that it appears to be a result of Iapetus' 
synchronous rotation about Saturn, which leads to un- 
equal radiation and micrometeorite bombardment fluxes 
on the leading and trailing hemispheres of the satellite. 
This special geometric circumstance is presumably not 
relevant to the case of small bodies in heliocentric orbit. 

In the outer solar system, lower temperatures and 
greater distances make the detection of thermal radia- 
tion increasingly challenging, even w ith the most sens i- 
tive infrared satellites in space (e.g. iCruikshankl [2005) . 
Consequently, only the reflected lightcurve is available, 
and the interpretation must be based on the assumption 
that the surface albedo variation is minimal. As we will 
see, support for this assumption comes not only from the 
analogy with the (generally uniform) main-belt asteroids, 
but from the remarkably symmetric lightcurves displayed 
by most outer solar system objects. Rotational symme- 
try is expected for figures of equilibrium having uniform 
surface albedos but is not a natural consequence of sur- 
face albedo markings. 

The lightcurves of several large Kuiper Belt Objects 
(KBO s), notably (2000 0) Varuna (I Jewitt fc Sheppardl 
[2002] ) and 2003 EL61 (jRabinowitz et al.l 120061 ). sug- 
gest that these are high-angular momentum bodies 
in which the shape has been deformed by rapid ro- 
tation. Other objects may be contact binary sys- 
tems, as has lo ng been suggested for Jovian Tro- 
jan (624) Hektor (lHartmann fe Cruikshanld 119781 119801 : 
IWeidenschilHnlll980j) and recently, for KBO 200 1 QG 298 
(jSheppard fc Jewittjf200l iTakahashi fc Ip|l2004D . These 
systems are interesting since, under conditions of rota- 
tional equilibrium, the period and the shape (both of 
which can be inferred from lightcurve data) are uniquely 
related to the bulk density. Lightcurves of these objects 
may thus be interpreted in terms of a fundamental geo- 
physical property that is otherwise difficult to measure. 

In this paper, we discuss the lightcurves of specific solar 
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TABLE 1 
(624) Hektor lightcurve data. 



JD a 


e b 




Am data d 






[dog] 


[dog] 


[mag] 


[mag] 


JD2435989 


74.9 


+4.4 


0.775 


0.737 


JD2438795 


24.8 


+4.1 


0.113 


0.063 


JD2439556 


52.5 


-5.3 


0.398 


0.302 


JD2439977 


86.3 


+4.1 


1.055 


1.048 



Note. — Sources cited in Section |3~T1 
a Date of observation 10 Aspect angle c Phase 
angle d Data lightcurve range Best-fit model 
lightcurve range 



TABLE 2 
List of objects to fit. 



Object a Family 13 H [mag] c D c [km] d P [hr] c Am [mag] f p[kgm" 3 ] g 



Triaxial ellipsoids 



2003 EL 6 i 

Varuna 
2000 GNiTi 


KBO 
KBO 
KBO 


0.2 
3.2 
6.0 


1450 
900 
360 


3.9154 
6.3442 
8.329 


0.28 + 0.04 
0.42 + 0.02 
0.61 + 0.03 


2585^ 
992+^ 

1Q46+1380 


Roche binaries 


2001 QG 2 98 
2000 GNiTi 
(624) Hektor 


KBO 
KBO 
Trojan 


6.9 
6.0 
7.5 


240 
360 
180 


13.7744 
8.329 
6.9225 


1.14 + 0.04 
0.61 + 0.03 
1.2 h 


590± 4 4 7 3 
2480:% 



Note. — Sources cited in the text. KBO 2000 GN171 is intentionally listed twice, as its nature 
is uncertain. 

a Object designation. 13 Object family. c Absolute magnitude. d Approximate equivalent circular 
diameter. Rotation period.' Peak-to-peak lightcurve range. 8 Estimated density. 11 Maximum 
predicted amplitude, at 9 = 90 deg. 



system bodies in terms of rotational equilibrium models, 
paying particular attention to high angular momentum 
systems and contact binaries. Our models address the 
effects of the surface scattering on the derived system 
parameters. Prototype contact binary (624) Hektor is 
examined in detail, taking advantage of voluminous high 
quality data published for this object over a range of 
aspect angles (see Table [1]). The models are then applied 
to four well-observed KBOs (Table [2]) and used to place 
quantitative contraints on their properties in a consistent 
formalism. Indeed, the uniformity of approach is one of 
the strengths of our simulations. 

2. LIGHTCURVE SIMULATIONS 

2.1. Jacobi Ellipsoids 

The formalism associated with the ellip soidal figures of 
equilib rium is described in great detail in IChandrasekharl 
(1969). A homogeneous, fluid body spinning in free space 
will assume a shape that balances self-gravity and the 
inertial acceleration due to rotation. This means that 
the triaxial shape of such a body is a function of its 
spin frequency and density. The equilibrium figures of 
isolated, rotating bodies are the Maclaurin spheroids and 
the Jacobi ellipsoids. The former are oblate spheroids 
and the latter are triaxial ellipsoids, and in both cases 
the rotation is about the shortest physical axis. We are 
interested only in the Jacobi ellipsoids because oblate 
spheroids have rotational symmetry and thus produce 



flat lightcurves. 

The shapes of Jacobi ellipsoids in ter ms of the semi- 
axes (a, 6, c) can be obtained by solving (|Chandrasekharl 
1969) 



2 j 2 

a o 



1 



o (a 2 + u) (6 2 + u) A(a, 6, c) 

■2 



du 



(1) 



o (c 2 + u) A(a, b, c) 



du 



where A (a, b, c) = yj (a 2 + u) (b 2 + u) (c 2 + u). The spin 
frequency u> and density p are related to the shape by 

,2 

2abc I — ^77^- — -du (2) 



nG p 



(a 2 + u) (b 2 + u) 



where G is the gravitational constant. We solved Eq. |(3J) 
for values of b/a between 0.43 and 1.00 in steps of 0.01, 
and used the solutions, together with Eq. ([2]), to cal- 
culate lo 2 /(ttG p) which relates the spin period to the 
body density for each of the equilibrium triaxial ellip- 
soids. Figures with b/a < 0.43 are unstable to rotational 
fission (|Jeandll919f ). 

The derived shapes are then raytraced at regular inter- 
vals spanning a full rotation period. This produces a set 
of frames from which the lightcurve is extracted by inte- 
grating the total light in each one of them. In this way 
we generate a database of lightcurves of figures of equi- 
librium which can be used to compare to the lightcurve 
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-0.001 
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0.001 
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phase 

Fig. 1. — Accuracy of raytraced versus analytical solution, for 
an ellipsoid with a "lunar" type surface, an axis ratio b/a=2/3, 
observed at aspect angle 9 = 90° and phase angle a = 0° . Bottom 
points and right ordinate axis show difference between the points 
(raytraced) and the line (analytical) shown at the top. 

data. As described below we run our simulations for two 
surface scattering laws. 

The raytracing is done using the main en- 
ine of the open source software POV-Ray 
http://www.povray.org[ ). The surface scattering 
routines were rewritten to permit accurate control 
of the scattering function. To test the accuracy of 
our raytracing method we simulated the lightcurve of 
a triaxial ellipsoid with an axis ratio b/a, observed 
equator-on (8 = 0°) at zero phase angle (a = 0°), usin g 
a "lunar" surface scattering function (see Section I2.3() , 
and compared the result with the analytical solution for 
the same configuration, given by 



m(b/a, <j>) = 2.5 log 10 a/1 + [(b/a) 2 - 1] cos 2 (2 7r</>). 

(3) 

where 4> S [0, 1] is the rotational phase of the ellipsoid. 
The result is plotted in Fig. [TJ The raytraced lightcurve 
deviates ~ 0.1% from the analytical solution which is 
negligible when compared to the uncertainties in the pho- 
tometric data, typically ~ 2 to 3% on real astronomical 
objects. 

2.2. Roche Binaries 

Lightcurves from an eclipsing binary asteroid consist- 
in g of two spheres in circ u lar or bit have been presented 
by IWiiesinghe fc Tedescol (| 1979( 1. As the binary separa- 
tion becomes comparable to the scale of either compo- 
nent, mutual gravitational forces will deform the bodies, 
increasing the lightcurve range over the ma ximum (fac- 
tor o f two) possible for equal-sized spheres (|Leone et al.l 
1984). To model this mutual deformation in close and 
contact b inary systems we use the Roche binary ap prox- 
imation (jChandrasekharl 119631 : iLeone et all 1 1984( 1. In 
this approximation each component is considered to be 
a Roche ellipsoid, which is the equillibrium shape of a 
satellite orbiting a spherical, more massive primary. The 
tidally deformed shape of the secondary is assumed to be 
solely caused by the spherically symmetric gravitational 
gradient due to the primary. Each component's shape is 
calculated separately using reciprocal values, q and 1/q, 
for the mass ratio. Clearly, such an approximation in- 
troduces the most error when calculating shapes of close 
binaries with mass ratios near 9=1. In these situations 
the elongation of the binary components is underesti- 
mated which leads to smaller lightcurve ranges. With 



the further assumptions that the binary is tidally locked, 
and that the components have equal density p, and orbit 
the center mass in circular paths, the mass ratio q, the 
shape of one of the components (b/a, c/a) . and the orbital 
frequ ency u> can be calculated by solving (|Chandrasekharl 
19G3) 



(3 + 1/q) a 2 + c 2 _ a 2 A 1 - c 2 A 3 



(l/q)b 2 

, ,2 



1 + q -nGp 
A x , A 2 , 



2abc- 



b 2 A 2 
a 2 A 1 



c 2 A 3 
c 2 A 3 



(3 + 1/q) a 2 

with A\ , A^ and A 3 

(|Chan drasck har fc Lebovit3ll962fl 

M = . 2 , , — L [F(0, <f) - E(9, <f>] 



given 



a 3 sin 3 ( 



sin 



A 2 = 



1 



a 3 sin 3 ( 



sin 9 cos 2 9 



■x 



E(6,<t>)-F(d,< 



A, 



cos z e 
b 



a 3 sin 3 i 



E(6, 



(4a) 

(4b) 
by 

(5a) 
(5b) 

(5c) 



where E(9, <f>) and F(9, 4>) are the standard elliptic inte- 
grals of the two kinds with arguments 



= arcsm 



b 2 



and 6 — arccos 



(!)■ < 6 > 



a^ — \a/ 
Each root of Eq. (|4"aj) . corresponding to a Roche ellip- 
soid solution, can be calculated by setting 3 of the 4 
parameters q and (a, b, c), and solving for the 4th by 
interpolation. For the primary, we set a = 1 and cal- 
culate b for each combination of q — q m i a , . . . , 1.00 and 
c = Cmin, ■ ■ • , 0.99, both in steps of 0.01. The procedure 
is repeated using 1 /q instead of q to calculate the shape 
of the secondary, (a',b',c'). A valid Roche binary solu- 
tion is obtained if two sets, (q, a, b, c) and (1/q, a' ', 6', c'), 
yield the same value uj 2 (7r G p)^ 1 when replaced into the 
right-hand side of Eq. (|4b[) . Table [3] shows the solutions 
for q = 0.25. 

2.3. Surface Scattering 

The surface scattering properties of KBOs are un- 
known. For this reason the amount of sunlight reflected 
from a KBO is usually taken to be proportional to its 
geometrical cross-section. However, the total range of 
the lightcurve of a convex object increases significantly 
if there is limb darkening. The two simplest scatter- 
ing laws generally used to model planetary surfaces are 
the Lommel-Seeliger and Lambert laws. The Lommel- 
Seeliger law is a single scattering model, suitable for low 
albedo surfaces, and can be conside red a simplification 
of the well known Hapke model (Li 2005|). The Hapke 
model is inappropriate for this work because it has many 
parameters which cannot realistically be constrained us- 
ing lightcurve data. To model relative brightness varia- 
tions, which is what is needed to generate lightcurves, the 
Lommel-Seeliger law requires no parameters: it depends 
solely on the cosines of the incidence and emission angles 
(i and e, the angles between the surface normal and the 
directions to the light source and to the observer). The 
Lommel-Seeliger reflectance function is thus 



fi + p a 



(7) 
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TABLE 3 

Roche binary solutions for mass ratio q = 0.25. 



q a 


B/A b 


C/A h 


b/a c 


c/a c 




l c 


r\ ornnn 

(J.zoUUU 


U.91o74 




0.51426 


r\ a or\r\r\ 

(J.480UU 


U.luozo 


1.19ZZZ 


U.zoUUU 


n noono 


U.o4UUU 


0.61175 


U.O f UUU 


U.lUlo/ 


l.Zol9o 


(J.25UUU 


n Aoom 

0.92891 


U.85UUU 


U.dd4v>o 


r» c nnnn 

U.oilKXJ 


U.(J9oo4 


1. ,5447a 


U-ZOUUU 




U.oDUUU 


U. I UOoZ 


U.OOUUU 






0.25000 


0.94037 


0.87000 


0.73534 


0.69000 


0.08368 


1.45810 


0.25000 


0.94584 


0.88000 


0.76469 


0.72000 


0.07755 


1.51789 


0.25000 


0.95114 


0.89000 


0.79333 


0.75000 


0.07109 


1.58466 


0.25000 


0.95629 


0.90000 


0.81200 


0.77000 


0.06558 


1.64436 


0.25000 


0.96128 


0.91000 


0.83935 


0.80000 


0.05864 


1.72902 


0.25000 


0.96612 


0.92000 


0.85713 


0.82000 


0.05282 


1.80713 


0.25000 


0.97081 


0.93000 


0.88311 


0.85000 


0.04550 


1.92216 


0.25000 


0.97537 


0.94000 


0.89996 


0.87000 


0.03943 


2.03365 


0.25000 


0.97979 


0.95000 


0.91642 


0.89000 


0.03328 


2.17018 


0.25000 


0.98407 


0.96000 


0.93250 


0.91000 


0.02706 


2.34450 


0.25000 


0.98824 


0.97000 


0.95588 


0.94000 


0.01924 


2.65422 



a Mass ratio b Primary 
squared in units of ir G 
tance in units of A + a 



axes ratios c Secondary axes ratios d Orbital frequency 
p, where G is the gravitational constant Orbital dis- 




0.4 0.6 0.8 
phase 

Fig. 2. — Lightcurves of a Roche contact binary at different phase 
angles and for different scattering laws. In the "uniform" case every 
point of the surface that is illuminated by sunlight reflects exactly 
the same amount of light back to the observer. Both primary and 
secondary have axis ratios b/a = 0.67 and c/a = 0.60, and the 
components are in contact. 

where /zo = cosi and /i = cose. 

The Lambert scattering law is a simple description of 
a perfectly diffuse surface. It assumes that a light ray 
that enters the material is multiply scattered and thus 
leaves the surface in a random direction. As such, it is a 
multiple scattering law which adequately describes high 
albedo surfaces. The Lambert reflectance function is 

r-L °c Mo- (8) 

The Lommel-Seeliger and Lambert scattering functions 
are taken here as representative of low albedo "lunar" - 
type surfaces and high albedo "icy" -type surfaces, re- 
spectively. 

Figures and [3] compare the lightcurves of Roche con- 
tact binaries using both lunar and icy scattering mod- 
els, at four different phase angles. A "uniform" scat- 
tering law is also plotted for comparison. The uniform 
model assigns equal brightness to any illuminated point 
on the surface, and thus represents the illuminated cross- 
section. Figure is for a contact binary with equal size 
components, while Fig. [3] represents the case of differ- 
ent sizes for primary and secondary (see Fig. [3]). A few 
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Fig. 3. — Same as Fig. [2] but for a binary with different size 
components. Mass ratio is q = 0.67 and the components have axis 
ratios (B/A = 0.77, C/A = 0.69) and (b/a = 0.53, c/a = 0.49). 




Fig. 4. — Three-dimensional rendering of the binary used to 
produce Fig. |3](a = 60deg). Rotational phase (0 = 0, 72, 144, 216, 
and 288 deg) runs from left to right and top to bottom rows show 
lunar, icy, and uniform surface scattering functions. 

conclusions can immediately be drawn from inspection 
of Figs. Hand O 

• At low phase angle the lunar model produces neg- 
ligible limb darkening and is thus equivalent to the 
simpler uniform model. Only at large phase angles 
(> 30 deg) does the uniform approximation fail to 
follow the lunar scattering law. Icy surfaces, how- 
ever, always produce larger lightcurve ranges for 
the same shape, which implies that assuming uni- 
form scattering when interpreting lightcurves of icy 
objects will tend to exaggerate the inferred shape 
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elongation. 

• The lightcurve minima become broader with in- 
creasing phase angle. This fact implies that V- 
shaped minima are strictly diagnostic of a close 
binary configuration only when viewed near zero 
phase angle. Fortunately, this is necessarily the 
case for all observations of Kuiper belt objects. 

• Observations at different phase angles shift the 
minima and maxima of the lightcurve in rotational 
phase. This effect should be taken into account 
when fitting a single spin period to observations 
taken at different phase angles. For instance, fail- 
ure to fit all the data with a single period does 
not necessarily imply complex (non-principal axis) 
rotation when observations over a wide range of 
phase angles are compared. 

In our simulations we use the lunar- and icy-type sur- 
face laws separately to assess how different surface prop- 
erties affect our results. As shown below, we find that 
different surface properties do not significantly change 
our density estimates. It is also shown that while in some 
cases the choice of a particular scattering law clearly im- 
proves the lightcurve fit, in other cases the result is de- 
generate as far as surface properties are concerned. The 
simulations further assume that the surface albedo is uni- 
form and the same for both components of the binary. 

2.4. Observational Geometry 

Owing to their large distances from the Sun and the 
Earth, KBOs can only be observed at small phase an- 
gles (a < 2deg). (20000) Varuna is known to exhibit a 
pronounced oppos ition effect at phase angle a < 0.1 deg 
(|Hicks et al.ll2005T ) w hich seems to affect t he extent of its 
brightness variation (jBelskava et al.l [ 2006| ). In the range 
0.1 deg < a < 2.0 deg the phase cur ve of (20000) Varuna 
is linear (jSheppard fc Jewittl 120021 ) and the lightcurve 
unnaffected by phase effects. We place our simulations 
well within this linear regime by using a phase angle 
a = 1 deg. In addition, since most of the data being 
modeled have a > 0.1 deg the opposition effect is likely 
unimportant to the conclusions presented here. 

The aspect angle is the angle between the line-of-sight 
and the spin axis of the system. We simulated lightcurves 
for two values of the aspect angle: 90 degrees (equator- 
on) and 75 degrees. Given that the spin axis orientations 
are unknown, the models with a tilted orientation allow 
us to investigate how the aspect angle affects our conclu- 
sions. In the case of (6 24) Hektor the spin (orb ital) axis 
orientation is known (|Dunlap fc Gehreisl Il969f ) and we 
were able to perform simulations using the actual aspect 
and phase angles at the moment the data were taken (see 
Table [[]). 

3. BEST FIT SOLUTIONS 

3.1. (624) Hektor 

Jovian Trojan (624) H ektor has long been recognized as 
a likely contact binary (lHartmann fc Cruikshanklll978l 
119801 : rWeidenschillinglll980T ). Numerous lightcurve data 
spanning a long time base (1957-1968) have been col- 
lected for this object, which allowed the determina- 
tion of the pole orientation (|Dunlap fe Gehreisl 119691 : 
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Fig. 5. — Quality of fit as a function of scattering function (top- 
left), mass ratio of binary components (top-right), orbital distance 
(bottom-left), and bulk density of binary components (bottom- 
right) for (624) Hektor. To avoid cluttering, only Roche models 
with mass ratio values, q, multiple of 0.05 are plotted. 
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JD2435989 V ' 
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Fig. 6. — Roche binary model lightcurve superimposed on data 
for (624) Hektor taken at 4 different aspects. 



Ide Angelii Il995f ). Depending on the orbital configu- 
ration, the lightcurve range of (624) Hektor varies be- 
tw een 0.1 mag and 1.2 m ag (see T able [B ■ As note d 
by IWeidenschilThii (|1980l see also iLeorie et all Il984h . 
lightcurve ranges above 0.9 mag cannot be produced by 
rotation of a single equilibrium figure and, instead, are 
most likely produced by a tidally distorted contact bi- 
nary. Besides the large range of variability, the lightcurve 
morphology of (624) Hektor exhibits V-shaped minima 
and round maxima, also characteristic of tidally de- 
formed, contact binary systems. 

Making use of the ex tensive dataset presented in 
iDunlap fc Gelireli (fl969l) . we have selected the Roche 
binary model (see Section |2~2|) that simultaneously best 
fits the observations at four observing campaigns (see Ta- 
ble[T|). The quality of fit is measured by xVXbesf i- e -' ^ ne 
ratio of the chi-squared value of each model to the chi- 
squared value of the best fit model. This corresponds to 
a reduced chi-squared, \ 
fit model has Xrcd = !■ 
associated with the data for (624) Hektor are not known 
with certainty, which does not allow us to reliably calcu- 
late the reduced chi-squared for each model. In Table [4] 



2 , d , if one assumes that the best 
We do this because the errors 
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TABLE 4 

Best three Roche binary model fits to (624) Hektor. 



SL a 


Q h 


B/A c 


C/A c 


b/a d 


c/a d 


u, 2 /(nGp)° 


d/(A + a) f 




Y 2 /Y 2 h 
A / Abcst 


Lunar 


0.62 


0.80 


0.72 


0.47 


0.43 


0.122 


0.98 


2480 


1.00 


Lunar 


0.67 


0.77 


0.69 


0.53 


0.49 


0.128 


1.00 


2374 


1.06 


Lunar 


0.65 


0.79 
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Fig. 7. — Rendering of best-fit Roche binary model for (624) Hektor for the 4 geometries listed in Table [T] Rotational phase (<j> = 
15, 60, 105, 165, 210, 255, 315, and 360 deg) runs from left to right while aspect angle decreases from top to bottom. 



we show the three Roche models that best approximate 
(624) Hektor's lightcurve. Figure [5] shows how well our 
simulations are able to determine the density, orbital dis- 
tance, mass ratio, and surface scattering properties of the 
(624) Hektor system. In the upper left panel we plot his- 
tograms of quality of fit for the two scattering laws con- 
sidered. A lunar scattering law clearly produces better 
fits, which is no surpr ise as (624) Hektor is known to have 
a low albedo (~0.06, ICruikshan"klll977t iFernandez et all 
2003). In the lower left panel the orbital distance of the 
Roche binary models is plotted versus quality of fit. The 
orbital distance is given in units of the sum of the semi- 
major axes of the primary, A, and of the secondary, a. 
The minimum is centered around d/(A + a) — 1, which 
corresponds to the binary components being in contact. 
Values d/{A + a) < 1 arc unphysical but we have de- 
cided to keep them for two reasons. Firstly, because they 
may result from th e approximations of the Roche model 
(jLeone et al.lll984l ). The orbital distance is calculated us- 
ing Kepler's 3rd law assuming point masses for the binary 
components. Since these have elongated shapes, gravity 
will be enhanced meaning the Roche model may underes- 
timate the orbital distance. Secondly, because the mech- 
anism that brought the two components together and 
formed the binary may have produced some deformation 
(a "crush" zone) around the point of contact, bringing 
the objects closer together than the situation of contact 
between two perfectly hard ellipsoids. Considering mod- 
els with xVXbest < 2, which roughly corresponds to a 
1-cr criterion, we find that the two components of (624) 
Hektor are separated by d/(A + a) = LOO^Sgg, which 
is consistent with contact. The top right panel of Fig. [5] 
shows that the model poorly constrains the mass ratio of 



the binary. The mass ratio which corresponds to the best 
fit Roche binary is q = 0.62, but the range of q values that 
fall within 1-cr of the best fit is broad. As for the bulk 
density of (624) Hektor (bottom right panel of Fig. [5]), 
perhaps the most interesting quantity to come out of our 
simulations, we find p = 2480^292 kgm -3 . This value 
closely confirms an ea rlier estimate: p ~ 2500 kgm -3 
(|Weidenschillinglll980f ). 

Figure [6] sh ows how the best Roche b inary model com- 
pares to the iDunlap fc Gehrela (|1969t , lightcurve data. 
The differences in lightcurve range from one campaign 
to the next reflect the effect of the observational geom- 
etry. Given the simplicity of the model, the agreement 
is remarkable and lends strong support to the idea that 
(624) Hektor is a contact binary. The model for Ju- 
lian day 2438795 shows the largest departure from the 
data. This is to be expected given the small aspect an- 
gle, 9 = 24.8 deg. The cross-section of the binary varies 
little at such unfavorable geometry and brightness vari- 
ations must be attributed to irregularities on the surface 
of the Trojan, which are not accounted for in the model. 
In Fig. [7] we show the best fit model rendered at the four 
aspect angles and for eight values of rotational phase. 
Recent observations using the LGS AO system at the 
Keck-II telescope suggest that (624) H ektor may have a 
bilobated shape (|Marchis et al.ll2~006bf ) and lend further 
support to the results presented here. 

3.2. 2001 QG 29 8 

Kuiper belt object 2001 QG298 completes a full rota- 
tion every P = 13.77 hr and its brightness var ies by 
Am = 1.14 ± 0.04 mag (jSheppard fc Jewittl l200l . The 
large range of brightness variation and relatively slow ro- 
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TABLE 5 

2001QG298 MODEL FIT. 



S.F. a 9 h q c B/A d C/A d b/a c c/a c u 2 / (n G p) f d/(A + a)S p h X 2 /xi 



Jacobi ellipsoid 



Lunar 


90 




0.43 


0.34 






0.283 




271 


2.21 


Icy 


90 




0.56 


0.41 






0.327 




234 


2.50 


Lunar 


75 




0.43 


0.34 






0.283 




271 


6.60 


Icy 


75 




0.50 


0.38 






0.310 




248 


2.51 


Roche binary 


Lunar 


90 


0.84 


0.72 


0.65 


0.45 


0.41 


0.130 


0.90 


590 


1.00 


Icy 


90 


0.44 


0.85 


0.77 


0.53 


0.49 


0.116 


1.09 


659 


1.09 


Lunar 


75 


1.00 


0.44 


0.40 


0.44 


0.40 


0.135 


0.76 


568 


1.62 


Icy 


75 


0.73 


0.74 


0.67 


0.54 


0.49 


0.130 


0.98 


589 


1.16 



^ Surface type. Aspect angle. C Mass ratio of the binary components. 

S h 

frequency of triaxial ellipsoid (binary). Binary orbital separation. 

of best-fit model. 
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Fig. 8. — Quality of fit versus bulk density (top-left) and axis ra- 
tios (top-right) of Jacobi ellipsoid models, and versus bulk density 
(bottom-left) and orbital distance (bottom-right) of Roche binary 
models for 2001 QG298 lightcurve data. Jacobi ellipsoid models are 
plotted for all four combinations of surface properties and observa- 
tional geometry listed in Table [5] while Roche binary models are 
plotted for both lunar- and icy-type surfaces at an aspect angle 
= 90°. 
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Fig. 9.— Models that best fit the lightcurve of 2001 QG298 for 
each combination of scattering law and geometry. Solid and dotted 
lines indicate lunar and icy surface scattering models, respectively. 
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Fig. 10. — Rendering of best-fit Roche binary 

model for 2001 QG298 . Rotational phase (<f> = 

45,90,135,165,210,240,285,315, and 360 deg) runs from left 
to right and top to bottom. 

tation provide compelling evi dence that 2001 QG298 is a 
contact or very close binary (|Sheppard fc JewittJ 12004'). 
We attempted to fit both Jacobi ellipsoid and Roche bi- 
nary models to the lightcurve data on 2001 QG298, using 
the two surface scattering laws and the two different ob- 
servational geometries described in Sections 12.31 and 12.41 
Table [S] and Figures [5] and [5] present a summary of the 
best fit models for different combinations of scattering 
law and observational geometry. 

Clearly, the Roche binary simulations fit the data bet- 
ter. Furthermore, the binary model favored by the data 
has a lunar-type surface and an equator-on geometry (see 
Fig. [9j. However, choosing an icy- type surface does not 
result in a significantly poorer fit. Indeed, the lower left 
panel of Fig. [9] suggests that an intermediate scatter- 
ing law is needed to fit the shallower minimum in the 
lightcurve data. Models tilted 15 degrees toward the line 
of sight are unable to fit the deeper V-shaped minimum 
in the data. Taking all the Roche simulations into ac- 
count we find that 2001 QG298 should have an orbital 
separation d/(A + a) — 0.90^oi4 (contact binary) and 
a bulk density p = 59O+47 3 kgm~ 3 . The uncertainty in- 
tervals are established in the same way as was done for 
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Fig. 11.— Same as Fig. |U but for 2000 GN m lightcurve data. 
Jacobi ellipsoid models are plotted for all four combinations of 
surface properties and observational geometry listed in Table [6] 
while Roche binary models are plotted for lunar-type surface and 
aspect angle 9 = 75°. 




Fig. 12.— Same as Fig.EJbut for 2000 GNm ■ 



(624) Hektor (see Section l3Tj) . Inspection of Tableland 
Fig. [5] shows that the best icy-type surface models have 
d ~ 1.09 and p ~ 660 kg m -3 . Given that the chosen 
surface scattering laws represent extreme cases (the sur- 
face of 2001 QG 2 98 probably combines single and multiple 
scattering behaviour) , we must conclude that the density 
we find does not depend strongly on the specific surface 
scattering properties of the KBO. The same applies to 
the binary components being in (or very close to) con- 
tact. The best fit Roche binary is rendered in Fig. [T0l 
Our results are in good agreement with an independent 
but similar attempt t o fit this object's ligh tcurve with a 
Roche binary model (|Takahashi fe Ip||200l ). 

3.3. 2000 GN m 

The rotational properties of KBO 2000 GN m 
also make it a good candidate contact binary 
(jSheppard fc Jewittll2004h . Its spin period and lightcurve 
range are P = 8.329 hr and Am = 0.61 ± 0.03 mag 
(|Sheppard fc Jewittll2002l) . However, as can be seen in 
Table [B] and Figures [TT] and the lightcurve fitting re- 
sults are not definitive about the nature of this KBO. A 
Roche binary solution is the one that best fits the data 



Fig. 13.— Same as Fig.[l0]but for 2000 GN m . 
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Fig. 14.— Same as Fig. |8] but for 2003 EL 6i lightcurve data. 
Jacobi ellipsoid models are plotted for all four combinations of 
surface properties and observational geometry listed in Table [7] 
while Roche binary models are plotted for icy-type surface at an 
aspect angle 9 = 75°. 

(see Tabled but it is not significantly better than a single 
Jacobi ellipsoid model. Inspection of the lightcurve fits 
(Fig. H"2"]) suggests that while the Jacobi ellipsoid model 
follows better the overall shape of the lightcurve, it is 
not able to reproduce the different minima present in 
the data. A Roche binary solution produces different 
minima, but seems to require an aspect angle 6 < 90 deg 
and a low mass ratio (q ~ 0.3) to be able to reproduce 
a lightcurve range as low as Am = 0.61 mag. The pre- 
dicted orbital separation is d/(A+a) = 1.09lto io* Lunar- 
type surface models consistently produce better fits than 
icy-type, irrespective of the nature (Jacobi ellipsoid or 
Roche binary) and orientation of 2OOOGN171. The in- 
ferred density is model dependent, but for each of the 
two models it does not depend much on scattering prop- 
erties nor on geometry. If 2000 GN 171 is taken to be a 
binary, its density should be p ~ 2000 kg m -3 . Instead, 
if it is an elongated ellipsoid it should have a bulk den- 
sity p ~ 650 kg m -3 . The best Roche binary model for 
2000 GN m is rendered in Fig. [TJ 

3.4. 2003 EL 61 

The lightcurve of 2003 ELgi indicates a rotation pe- 
riod of P = 3.9 hr and a lightcurve tota l range of 
Am = 0.28 ± 0.04 (|Rabinowitz et al.ll2006f) . The ex- 
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TABLE 6 

2000GNl71 MODEL FIT. 



S.F. 



B/A d C/A d b/a c c/a c ^/(irGp) 1 d/{A + a)> 



Jacobi ellipsoid 



Lunar 


90 




0.62 


0.44 




0.342 




613 


1.05 


Icy 


90 




0.75 


0.50 




0.362 




579 


1.19 


Lunar 


75 




0.55 


0.41 




0.325 




645 


1.04 


Icy 


75 




0.71 


0.48 




0.357 




587 


1.20 


Roche binary 


Lunar 


90 


0.25 


0.92 


0.83 


0.51 0.48 


0.106 


1.19 


1972 


1.09 


Icy 


90 


0.25 


0.94 


0.87 


0.74 0.69 


0.084 


1.46 


2504 


1.66 



Lunar 
Icy 



75 
75 



0.34 
0.25 



0.89 
0.92 



0.81 
0.84 



0.45 
0.61 



0.42 
0.57 



0.108 
0.101 



1.09 
1.28 



1946 
2067 



1.00 
1.05 



ratios of primary. Axis rati< 



ndary' Spin (orbital) 



frequency of triaxial ellipsoid (binary). Binary orbital scparatn 
of best-fit model. 



3ulk density of the bodies (in kgm" J ). Ratio of X of model to X 



TABLE 7 
2003 ELgi model fit. 



S.F. a 






B/A d 


C/A d b/a c c/a e 




d/{A + a)s 




v 2 /v 2 1 




Jacobi ellipsoid 


Lunar 


90 




0.80 


0.52 


0.367 




2585 


1.00 


Icy 


90 




0.88 


0.55 


0.372 




2551 


1.01 


Lunar 


75 




0.76 


0.50 


0.363 




2611 


1.02 


Icy 


75 




0.86 


0.54 


0.371 




2557 


1.00 


Roche binary 


Lunar 


90 


0.25 


0.98 


0.94 0.90 0.87 


0.039 


2.03 


24049 


6.08 


Icy 


90 


0.25 


0.99 


0.97 0.96 0.94 


0.019 


2.65 


49286 


6.45 


Lunar 


75 


0.25 


0.95 


0.89 0.79 0.75 


0.071 


1.58 


13339 


2.56 


Icy 


75 


0.25 


0.97 


0.93 0.88 0.85 


0.045 


1.92 


20841 


2.46 



Surface type. Aspect angle. Mass ratio of the binary comp. 
frequency of triaxial ellipsoid (binary). 8 Binary orbital separati. 
best-fit model. 



!. Axis ratios of primary. Axis ratios of the secondary Spin (orbital) 
Bulk density of the bodies (in kgm~ 3 Ratio of X 2 of model to X 2 of 
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Fig. 15.— Same as Fig. |9]but for 2003 ELe 



tremely fast rotation of 2003 ELgi implies that it must 
have a high density. Using a hydrostatic equilibrium 




Varuna 




Fig. 16. — Side (top) and tip (bottom) views of the Jacobi 
ellipsoid models of 2003 EL 6i (left) and (20000) Varuna (right). 
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Fig. 17.— Same as Fig. |U but for (20000) Varuna lightcurvc 
data. Jacobi ellipsoid models are plotted for all four combinations 
of surface properties and observational geometry listed in Table [8] 
while Roche binary models are plotted for lunar-type surface and 
aspect angle 6 = 75°. 

criterion, Rabinow itz et all (|2006f ) estimate p ~ 2600 — 
3340 kg m -3 . A binary solution would require a consider- 
ably higher (and unrealistic) density than a rotationally 
deformed ellipsoid. Binarity is also unlikely given the 
small range of brightness variation: for a binary to pro- 
duce such a shallow lightcurve, the pole axis must nearly 
coincide with the line of sight. Indeed, we find that no 
Roche binary model is able to satisfactorily fit this ob- 
ject's lightcurve data (see Table [7] and Fig. [14]). In the 
case of Jacobi ellipsoid models, all possible combinations 
of surface properties or orientation fit the data equally 
well. This is partly due to the large scatter present 
in the lightcurve data. However, the predicted density 
(p = 2585^44 kg m -3 ) depends little on specific choices 
of surface and geometry, and i s consistent with the p = 
2600 - 3340 kg m" 3 estimate of iRabinowitz et al.l (|2006f) . 
We find that the axis ratios of 2003 ELgi should fall in 
the ranges b/a = 0.76 - 0.88 and c/a = 0.50 - 0.55. The 
icy scat tering law (with 9 = 75 d eg; see Table [7J) is prefer- 
able, as IRabinowitz et all (|2006l ) find that 2001 EL 6 i has 
a high albedo (>0.6). The best Jacobi ellipsoid represen- 
tation of 2001 EL6i is shown in Fig. [16] 

3.5. (20000) Varuna 

The rotational properties of this object (P = 6.34 hr 
and Am = 0.42 ± 0.02) were interpreted in the con- 
text of ellipsoidal figures of e quilibrium and a density o f 
p ~ 1000 kg m" 3 was derived (| Jewitt & Sheppardll2003 ). 
Our simulations lend support to this result by showing 
that (20000) Varuna's lightcurve is well fit by a Jacobi el- 
lipsoid model. This is apparent from Table [H] and Figs.fTTl 
and |T5] Figure [TH] depicts the Jacobi ellipsoid model of 
Varuna. As in the case of 2003 ELgi, the quality of fit 
is degenerate as far as surface properties and orientation 
are concerned. Thus, depending on particular choices of 
these properties, Varuna's axis ratios lie in the ranges 
b/a = 0.63 - 0.80, c/a = 0.45 - 0.52. The bulk den- 
sity determination is again much more robust: we find 
p = 992±f^ kgm- 3 . 

4. DISCUSSION 

The photometric lightcurve of 2003 ELgi exhibits 
asymmetries which are not reproduced by the simple Ja- 
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Fig. 18.— Same as Fig.[9]but for (20000) Varuna. 



cobi ellipsoid model. Given the large size of this object 
(D cq ~ 1450, see Table [2]) which safely puts it in the 
gravity regime, we do not expect such irregularities in 
the lightcurve to be due to an irregular shape. Instead, if 
the lightcurve features are real, they could have the same 
origin as Pluto's brightness variation: albedo patches 
across the object's surface. Like Pluto, 2003 ELgi is large 
enough to hold a thin atmosphere, which might condense 
on the surface and cause the patches. 
The high density derived for (624) Hektor stands 



in contrast to the low value (p 



800 



+200 
-100 



kg m 



-3\ 



derived for resolved Trojan binary (617) Patroclus 
(jMarchis et al.l l2006a[) . The sizes of these two Trojans 
are similar: Hektor is 102±2 km in radius while Patro- 
clus is 70±2 km in radiu s, when measured and inter- 
preted in the same way (|Fernandez et al.l I2003D . The 
low density of Patroclus requires substantial porosity 
and also suggests an ice-rich composition. Hektor's 
density is consistent with zero porosity and a smaller 
or negligible ice fraction. This difference is puzzling, 
given that the albedos (0.057±0.004 and 0.050±0.005, 
respectively) are very similar, as are the optical reflec- 
tivity gradients [11.6±1 %/10Q0 A ( [Sawyer| I1991D and 
8. 8±1%/1000A dJewitt k Luulll990l ). respectively]. 

iMarchis et alJ (|2006al ) argued that the low density and 
inferred porous, ice-rich composition of (617) Patroclus 
was an indication that it originated in the outer part 
of the solar system. The high density of (624) Hektor is 
hard to explain in this context; does it indicate that (624) 
Hektor did not form in the outer solar system? Analo- 
gously, the similar size and surface properties of these 
two Jovian Trojans could be used to infer a common ori- 
gin. A similar density argument was used for the Satur- 
nian irregular satellite Phoebe (~ 220 km in radius), but 
in the opposite direction. The high density of Phoebe 
(p = 1630 ± 33kgm~ 3 ), when compared with that of 
other (regular) moons of Saturn, has been intepreted as 
indicative of an outer solar system origin on the basis 
that it matches the density of Pluto ([Johnson fc Luninel 
2005). As the examples above show, it is difficult to es- 
tablish a simple relation between formation region and 
bulk density - there may be no such relation - and there- 
fore the density of a body alone should not be used to 
infer its origin. 
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TABLE 8 
(20000) Varuna model fit. 



S.F. 



B/A d C/A d b/a c c/a c uj 2 /{nGp) f d/(A + a)£ 



Jacobi ellipsoid 



Lunar 


90 




0.69 


0.47 




0.354 




1020 


1.16 


Icy 


90 




0.80 


0.52 




0.367 




985 


1.00 


Lunar 


75 




0.64 


0.45 




0.346 




1045 


1.17 


Icy 


75 




0.77 


0.51 




0.364 
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1.00 


Roche binary 


Lunar 


90 


0.25 


0.92 


0.84 


0.61 0.57 


0.101 


1.28 
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5.48 


Icy 


90 
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Fig. 19. — Log density versus log equivalent circular diame- 
ter for the four KBOs modeled. Jacobi ellipsoid (Roche binary) 
fits are indicated by single (double) ellipsoid symbols. Pluto and 
Charon are plotted for comparison. KBO 2000 GN171 is plotted 
twice (connected by a dotted line). References in the text. 



Figure [19] shows the KBO densities from our simu- 
lations versus equivalent c ircular diameter; Pluto and 
Charon (jPerson et al.l 1200(1 ) are plotted for comparison. 
The sizes of 2001 QG 298 and 2001 GN m were calculated 
from their absolute magnitude assuming a 0.04 albedo, 
and the error bars exten d the albedo t o 0.10 . The size of 
(20000) Varuna is from Uewitt et all (120011, and that of 
2003 ELqi is calculated using its mass (|Rabinowitz et al.l 
2006) and the density derived here (see also Table [2]). A 
trend of increasing density with size is clear. Such re- 
lation may be caused by (1) a difference in composition 
(ice/rock ratio), with bigger objects having larger rock 
fractions, or (2) a trend in porosity, with larger objects 
being more compacted than their smaller counterparts, 
likely due to larger internal pressure. Although the latter 
effect is certainly present, it is unclear if it is the dom- 
inant cause for the trend. This size-density relation has 
been noted before (Jcwitt 20021) and seems to be present 
in di fferent populatio ns, e.g., KBOs and planetary satel- 
lites (jJewittJIin printj ). 



5. THE FRACTION OF CONTACT BINARIES IN THE 
KUIPER BELT 

Our simulations can be used to determine the 
lightcurve range of Roche binaries at arbitrary observ- 
ing geometries, and for two different surface types. We 
make use of this feature to refine an earlier estimate of 
the contact binary fracti o n in the Kuiper Be lt (KB; see 
ISheppard fc JewittJl2004f ). iLeone et aiT(|1984l ) considered 
that lightcurves with ranges between 0.9 and 1.2 mag- 
nitudes must be produced by tidally deformed contact 
binaries (see also IWeidenschiiling|[i980( ) . While the max- 
imum range of 1.2 mag is valid for lunar- type surfaces 
having negligible limb darkening, our simulations show 
that Roche binaries with icy-type surfaces (and thus sig- 
nificant limb darkening) can produce lightcurve ranges 
up to 1.57 mag. We searched our models for the binaries 
that produce these maximal lightcurve ranges (1.2 mag 
for lunar surface and 1.57 mag for icy surface) when ob- 
served equator-on, i.e., at aspect angle 9 = 90°. The 
aspect angle is measured between the line of sight and 
the pole axis of the binary. As the pole axis of such a 
binary moves away from the equator-on configuration (as 
the aspect angle approaches 0°) the lightcurve range be- 
comes smaller and smaller; let us denote by 6 min the as- 
pect angle at which the lightcurve range reaches 0.9 mag. 
If we assume that the pole axes of KB contact binaries 
are randomly oriented in space then the detected con- 
tact binary fraction is less than the true fraction by a 
(geometrical correction) factor cos6* m i n . 

Using our simulations we find the geometrical correc- 
tion factor to be ~ cos(81.4°) = 0.15 for lunar-type sur- 
faces and ~ cos(70.7°) = 0.33 for binaries with an icy- 
type surface (see Fig. . We use the fraction of 1 /34 
objects with large (> 0.9 mag) li ghtcurve range mea- 
sured bv ISheppard &: JewittJ (|2004[ ). as it constitutes the 
largest homogeneous survey for variability. Therefore, 
considering only lunar-type surfaces the true fraction of 
contact binaries is / ~ 1/(34 x 0.15) ~ 0.20. If we con- 
sider icy-type surfaces then we estimate the fraction to be 
/ ~ 1/(34x0.33) ~ 0.09. Two arguments make the latter 
of the two estimates a strong lower limit for the true frac- 
tion of contact (or close) KB binaries. Firstly, the two 
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Fig. 20. — Lightcurve range as function of aspect angle 6 for 
maximal Am Roche binaries with both lunar- and icy-type sur- 
faces. Top x-axis shows probability that the binary is observed at 
equal or larger 8. See text for details. 

surface types used here are simplified limiting models of 
how real planetary surfaces scatter light. Real objects 
presumably exhibit a degree of limb darkening between 
the two simulated here. Secondly, contact binaries with 
relatively low mass ratios produce shallower lightcurves 
which fall below the 0.9 mag threshold adopted here, 
and are not accounted for. Our estimate, new in that 
it includes the effect of surface scattering, substantiates 
the idea that a considerable population of contact /close 
binary objects in the Ku iper Belt may await discovery 
((Shepnard fc JewitiJ[20f)l . 

The Pan-STARRS all-sky survey 

( |http://pan-starrs.ifa.hawaii.edu/[ ) will scan the en- 
tire visible sky, down to mj;. ~ 24 on a weekly basis 
((Kaiser fc Pan-STARRS Teamll2005l ). Besides detecting 
all moving objects to that brightness limit, this cadence 
will allow (sparsely sampled) time series photometric 
studies, and thus the detection of high variability 
candidates, suitable for follow-up observations. The 
survey will therefore significantly improve the estimate 
of the contact binary fraction. The intrinsic fraction of 
KB contact binaries can provide impor tant constraints 
on b inary formation mechanisms (e.g. iGoldreich et al.l 
l2002f) and collisiona l evolution in the KB region 
((Petit fc Mousisll2004n . 

6. SUMMARY 

Mathematically unique interpretations of rotational 
lightcurve data are generally impossible. Nevertheless, 



lightcurves can, under physically plausible assumptions, 
convey invaluable information about the spins, shapes 
and densities of small solar system bodies. In this work 
we have explored the role of surface scattering proper- 
ties on the derivation of bulk densities from rotational 
lightcurves, using a quantitative model of rotationally 
deformed bodies. We find that: 

• With few exceptions, the choice of a particular scat- 
tering function does not strongly affect the densi- 
ties we obtain from our simulations. Instead, the 
presence of surface irregularities (lumps) and some 
albedo variegation (spots) on the objects sets the 
limit to the precision of our density estimates; sur- 
face lumps and spots make it impossible to find 
one idealized equilibrium shape that matches the 
lightcurve, leading to some degeneracy in the fits. 

• Our density estimates suggest a trend of increasing 
density with size. It is still unclear if such relation 
is mainly due to composition, to a trend in porosity, 
or to a combination of both. 

Confirming previous inferences, we find that: 

• The lightcurves of (20000) Varuna and 2003 EL 6 i 
are well-matched by rotational equilibrium mod- 
els in which the bodies are deformed by rota- 
tion into a triaxial shape. Jacobi ellipsoid models 
with uniform surface albedo and a range of limb 
darkening functions have been used to derive the 
bulk densities (Varuna: 992ti5kgm~ 3 ; 2003 EL 6 i: 
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• The lightcurves of Jovian Trojan (624) Hektor 
and KBO 2001 QG298 are well-described by con- 
tact binary models in which the densities are 
2480T.g u 2 kgm" 3 and 590t^ 3 kgm" 3 , respectively. 

• The high incidence of KBO lightcurves consistent 
with a contact binary interpretation suggests that 
these bodies are common in the Kuiper belt. 
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